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Abstract 



o 

Chaotic dynamics can be effectively studied by continuation from an anti-integrable 
limit. Using the Henon map as an example, we obtain a simple analytical bound on 
the domain of existence of the horseshoe that is equivalent to the well-known bound 
of Devaney and Nitecki. We also reformulate the popular method for finding periodic 
orbits introduced by Biham and Wenzel. Near an anti-integrable limit, we show that 
■ this method is guaranteed to converge. This formulation puts the choice of symbolic 

dynamics, required for the algorithm, on a firm foundation. 
AMS classification scheme numbers: 58F05, 58F03, 58C15 



1 Introduction 



The anti-integrable (AI) limit of a mapping is a singular limit in which the dynamics is 
not deterministic Q. Each orbit in the AI limit is given by a sequence of symbols, and the 
dynamics becomes the unrestricted shift operator on the symbols. It is well known that when 
the AI limit is "nondegenerate" many of the AI states can be continued away from the limit, 
becoming orbits of the original mapping [^, [|]. This continuation theory, a consequence of 
the implicit function theorem, is much simpler than the corresponding continuation from 
an "integrable limit" , which often requires the machinery of KAM theory. 

Biham and Wenzel 0] introduced a technique for finding periodic orbits of maps, in 
particular maps that are obtained from a variational principle. They generalized the notion 
of "gradient search" to find the minimum of the variational function by introducing a 
"pseudo- gradient" system of differential equations. This is obtained by multiplying the 
gradient by a diagonal matrix of signs. The basic idea is that this will allow one to find 
critical points other than the minimum. This method has become popular ||, |^, [8], || [l^. 
T5J 0, II, ||, in spite of the fact that it has no rigorous foundation. Indeed, 



Grassberger et. al [18] found examples where the method fails to find unique orbits for 



certain parameter values in the Henon map. 

'Useful conversations with Robert Easton, Bruce Peckham and Holger Dullin are gratefully acknowledged. 
JDM was supported in part by NSF grant number DMS-9623216. 
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In this paper we use the Henon map as an example to illustrate the usefulness of the 
AI limit. At the AI limit, the Henon map reduces to the full shift on two symbols, and 
we study orbits that continue from these states. We obtain an explicit bound on this 
continuation which happens to correspond exactly to the bound on the existence of the 
horseshoe obtained by Devaney and Nitecki [jl9| who used a geometrical argument. We 
then reformulate the Biham and Wenzel method and show that the pseudo-gradient system 
has a hyperbolic sink at the AI limit, and that this sink persists at least for the same 
parameter range for which we can continue the orbits. 



2 Anti-Integrable Limit 

The anti-integrable (AI) limit can be formulated in two ways. Consider a map, / : M — > M 
on a (I dimensional manifold M that depends upon some parameters. The general idea is 
to rewrite the system as an equivalent implicit relation F : M x M — > W d in such a way 
that F becomes singular when a parameter, say e, is zero For example if F[x,x') = 
eG(x,x') + H(x), then the "orbit" at e = corresponds to any sequence of zeros of if-the 
dynamics is not deterministic. We say that e = corresponds to the anti-integrable limit of 
the map /. If the derivative of H is nonsingular, then a straightforward implicit function 
argument can be used to show that (at least some of) the AI orbits can be continued for 
e ^ to orbits of the map / pC| , |5| . An AI limit with this property is called nondegenerate. 

Maps that are derived from a Lagrangian variational principle often have an AI limit. 
For example, orbits of the Henon map, written in the form 

x'\_fy-k + x 2 \ m 

y' J ~ \ -bx J ' [L) 

can be obtained as the critical point of an action, 

W[x] = Y,b- t L(x t ,x t+1 ), (2) 
t 

where we label time along the trajectory by a subscript t, and the factor 6~* allows the 
system to be non-area preserving. An orbit x = {. . . xt,x t +i, . . . } is a critical point of W, 
and yt = W^O^t-ii %t) = —Li(xt, xt+i). When the Lagrangian L can be put into the form 
L(x,x') = eT{x,x') — V(x), then e = corresponds to an AI limit ||. This form is quite 
natural for "mechanical systems" where T represents kinetic and V potential energy. 
For example, the Lagrangian for the Henon map is 

L(x,x') = —xx'-\ — x 3 — kx . (3) 
3 

Critical points of the action give the second order difference or Euler-Lagrange form of the 
map. To put (0) into a form appropriate for the AI limit, we rescale, defining 
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to obtain the new Lagrangian, 

z 3 

L{z,z') = e 3 L(x,x') = -ezz + — -z. (4) 

Critical points of the action correspond to sequences zt such that 

-e(z t+1 + bzt-i) +z?-l = 0. (5) 

Of course this implicit system is easily solved for zt+i whenever e ^ to give an explicit 
dynamical system. However, when e = (||) no longer defines a map, even though it still 
represents a critical point of W. This is the AI limit of the Henon map. At this limit, 
"orbits" correspond to any bi-infinite sequence z£§, where we denote the set of bi-infinite 
sequences with two symbols by 

S = {s: s t € {1,-1}, t £ Z}. 

Each s £ § continues to a unique orbit z(e) where z(0) = s when e is small enough 
|, §. In fact, we can obtain an explicit upper bound on the range of e for which this 
correspondence is guaranteed. It will be appropriate to use the l°° norm, HxHoo — sup^ 
and define Bm to be the closed ball of radius M around the point s, 

j Ba/(s) = {z:||z-s|| 00 <M}. (6) 

With this notation, we can prove: 

Theorem 1. For every symbol sequence s £ S, there exists a corresponding unique orbit, 
z(e), of the Henon map such that z(0) = s providing 



|e|(l + |6|) < 7oo = - 2/V~5 ~ 0.649839 . (7) 
The orbit z(e) is contained in the ball -Bmoo( s ) where 




7 7+ 7^ , (8) 

(9) 

To prove the theorem, we will write orbits z(e) of the Henon map as fixed points of an 
operator T whose t th component is 



T t (z) = sty/1 + e(z t+1 + bz t - 1 ), (10) 

where we choose the sign of the square root using a sequence s G S. When e = 0, the 
operator becomes T(z) = s, which trivially has a unique fixed point corresponding to the 
AI state, z = s. We will use the contraction mapping theorem to show that for small enough 
e the fixed point persists. The first step in the proof is to show that there is a domain C\ 
of the parameters (7, M) such that T is a contraction on Bm(s); this domain is illustrated 
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Figure 1: Contraction regions C n for 6=1 and n = 1,2, 30,oo. 

in Fig. |l|. This is improved by bounding the domain C n in the (7, M) plane where T n (z) is 
a contraction. It is easy to see that 

thus we can improve the bound on 7(e) by iteration of T. Figure [l| illustrates the parameter 
domains that we find for n = 1, 2, 30, and 00. The maximal 7 values for each of these iterates 
are denoted j n in Fig. ||. 

To demonstrate that T n is a contraction mapping we must show that for (^y, Tkf ) £ C n 

1. T" : B M (s) -> B M (s) 

2. I |Z?T^(z)| |oo < 1 for all z G Bju(s) 

The first requirement gives the solid curve in the figure defining a lower bound on M values 
in C n as well as right boundary where the range of T n is no longer real. The dashed 
boundary of C n is given by the second requirement. The maximal value of 7 is obtained in 
the limit n — >• 00. 

We begin with a lemma to prove the first requirement. 

Lemma 1.1. For any \j\ < 1/V2, and any M > ^00(7) there is an N such that T n : 
B M (s) -» B M (s) ,Vn>N. 

Proof: 

Let T be defined by (0). For any z G %(s) it is easy to see that 

a n <\\T n (m)\\ 0O <l3 n , (11) 
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where the sequences a n and (5 n are given by the iterations 



Pn+l = / (AO = y/l +lPn, 
OCn+X = \A _ lPn+1 , 

with the initial conditions (3q = 1 + M and oto = 1 — M. The map f(f3) has a single 
attracting fixed point 



7 + yV + 4 
Poo - 2 

Each of the a n must be real, so we must have 1 — 7/3 n > 0. This requirement gives a right 
boundary to the region in the (7, M) plane where T n exists. As n — > 00 these boundaries 
approach the vertical line defined by 



l-T&o = 0^7 = 1/V2, (12) 

which gives one of the bounds in the lemma. 
Finally, (11) implies that 

||T n -slloo < max(|a n - l\,\/3 n - 1|) = l-a n . 

Thus the requirement that T n maps Bm(s) into itself gives the implicit relation 1 — a n < M. 
As n approaches infinity, these domains approach the region 

M>Moo(7) = l-aoo. 

For each 7 the sequence M n for which l — a n — M n = converges monotonically to from 
above; therefore for any M > Moo (7), there is an N such that T n : %(s) — > Bm(s) , Vn > 
JV. □ 



Proof of theorem (|^) Let Bm(s) be defined by (g) and T by (p!o[). The map T n is a con- 
traction if I |DT n (z)| |oo < 1 for all z € Bm(s)- Using the chain rule gives 

||Z)T n ||oo< — f£ (13) 



From this and lemma 1.1, the map T n is a contraction in the region C n given by 

n 

C n = {(7, M) : 1 - a n < M , 7 < 2 «V" , Qj Real for j < n} . (14) 

i=i 

We find the boundary of by noting that the product in ( [T4] ) converges to aoo because 
the sequence otj converges geometrically to the fixed point. Thus there exists an N such 
that the map T n is a contraction for all n > N whenever 

7 < 2aoo • (15) 



3 THE METHOD OF BIHAM AND WENZEL 



6 



Using the form for we obtain the bound 

,2 



r < 4(i - 2/VE) , (16) 



which directly gives ([?]). This bound is clearly more restrictive than (12) so we can conclude 
that for any 

( 7 ,M) e Coo = {(7. AO : M > M oo{l), 7 < 7oo} 

there is an TV such that Vn > iV, the map T n (z) is a contraction map whenever z G Bm(s) 
The Banach fixed point theorem then implies that T n has a unique fixed point in Bm^ (s) 
for all n > N. But if z is a fixed point both of T^ +1 (z) and of T^ +2 (z) then it must be a 
fixed point of T. Since this is true for an arbitrary symbol sequence, for any 7 in the range 
dH), there is a corresponding unique orbit of the Henon map. □ 



By a completely different argument, Devaney and NiteckiJO]] proved that the non- 
wandering set of the Henon map is hyperbolic and conjugate to the 2-shift for exactly the 
same range of 7 that we give in Theorem |l|. 

3 The Method of Biham and Wenzel 

Biham and Wenzel [Q, |5| introduced a popular algorithm for finding periodic orbits of maps. 
They present the method in four steps 

1. Write the map in variational form with action W[x]. Critical points of W give the 
desired orbits. 

2. Use symmetries or some other technique to label the various critical points of W 
symbolically with a symbol sequence of signs, Sj. 

3. Choose an "arbitrary" initial guess for a periodic orbit, Xj(0),j = 0,1, ... n — 1 of 
period n. 

4. Introduce a pseudo-gradient dynamics with a fictitious time r to find an extremum of 
the action 

^ = -Diag(s)VlU[x] . (17) 

UT 

For a period n orbit, we set Xj +n (r) = xAt), and solve only the differential equations 
for x Q ,xi, . ..x n -\. 

Note that when Diag(s) = I this method is a standard gradient search for finding a local 
minimum of W , and if we set Diag(s) = —I, then we would simply be finding a local maxi- 
mum of W. More generally, the symbols are chosen in an attempt to "identify the critical 
points" of the action, but there is "no general way to find" the symbolic representation (BJ. 

The hope is that for each choice of Sj (and some range of map parameter) the system 
( |17| ) either has a "unique" attracting fixed point, or else the "corresponding" periodic orbit 
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does not exist. This claim seems at least plausible for the Henon map, since it has at 
most 2 n fixed points of period n, and there are precisely 2 n choices of periodic sequences s. 
However, Grassberger et. al [|l^] have found examples for the Henon map where the method 
fails to find unique orbits at certain parameter values. More generally, there is no guarantee 
that the fixed point of (|l7]) corresponding to the orbit of interest is attracting. Moreover, 
this system of equations certainly can have other attractors, including chaotic ones. 

A rigorous foundation for the Biham and Wenzel method can be given when the map 
has an AI limit. Suppose the discrete Lagrangian has the form L(x,x') = eT{x,x') — V(x), 
and that the potential V has a discrete set of nondegenerate critical points Crit(V) = 
{co, ci, C2, • • • Cm, c m+ i}, where we order the points, Cj < Cj + \ with the convention that 
Co = — oo and c m+ i = oo. When e = (|l7|) decouples, and becomes simply 

%—.,*(,,). 

Of course the equilibria of these differential equations are precisely the AI states, i.e., any 
sequence of critical points, Xj = Ck j , kj £ {1, 2, . . . m},j £ Z. Such an equilibrium is a sink 
if Sj = sgn(y"(cj- )). For example when the c/% are minima of V, then Sj = 1 as we expect. 
The basin of attraction of the sink is the box 

D s (0) = {x : c fe i < Xj < c fc +i} . 



Since hyperbolic fixed points of a vector field are preserved under a C 1 perturbation [21], 
we immediately have the theorem 

Theorem 2. If L = eT(x,x') — V(x) £ C 2 , and the critical points {ck} of V are nonde- 
generate, then there is an e max such that for all \e\ < e max and all sequences of critical 



points {ckj,j £ Z}, the system of differential equations (11), with Sj = sgn(V"(ckj)), has a 
hyperbolic sink, x*, that continues from Xj = c^. at e = . 



This reformulation of the Biham- Wenzel method points out several important things. 
First, the correct choice of the signs is determined by the signature of the potential V at the 
appropriate critical point. In fact, as originally formulated, the method cannot possibly find 
all orbits when V has more than two critical points. Rather, an orbit is determined first by 
the choice of anti-integrable state C&.; the sequence of signs Sj is determined subsequently. 



Finally, as we will explicitly demonstrate below, it is sensible to use the AI state as the 
initial condition for the method. 

For the Henon map, we use the scaled Lagrangian to obtain the system 

dz ■ 

-± = Sj (e(z j+1 + bzj-x) - z) + 1) . (18) 
When e = 0, the differential equations decouple and give the simple equations 

for which there is a unique attracting fixed point at Zj = Sj with a basin of attraction 

D s (0) = {z : -1 < Zj < oo} . 
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We can use theorem (|l]) to obtain a bound on e for the persistence of the hyperbolic fixed 
points of (|i~8|). We know from this theorem that there is a fixed point z* for each s when 
7 < Too- The linearization of ( |l8| ) around this fixed point is 

f = zl -1,0 • (19) 

Aij = Si [e(5 i+ ij + bSi-ij) - 2z*5ij] . 

By the Gerschgorin circle theorem, the eigenvalues A& of A are contained within the union of 
the disks centered at the diagonal elements with radius given by the sum of the magnitudes 
of the off-diagonal elements [22]. Since the row sum is bounded by |e|(l + |6|) = 7, we obtain 

A fc €|J{A:|A + 2|*;||<7}. 
3 

According to theorem ([[]), z* G -Bmoo(s), where M m is given by (^]); therefore, 

Re(X) < -2(1 -Moo) +7. (20) 
This is negative precisely when 7 < 700. Thus we have proven 

Theorem 3. When |e|(l + \b\) < 700, then the orbit of the Henon map labeled by the AI 
symbol sequence s is a hyperbolic sink for the system 



The basin of attraction of the sink initially includes the AI state z = s. We will next 
show that this remains true providing e is small enough. Thus the proper initial condition 
to use for the Biham-Wenzel equations is the AI state. 

Theorem 4. The AI sequence s is in the basin of attraction for the fixed point z* of 
providing 

|e|(l + |6|) < 0.555668. 

Proof: Suppose z* is an orbit of the Henon map for some fixed 7 < 700- Using dlg|), the 
deviation £ = z — z* from this orbit satisfies the system of equations 

f »4- (21) 

where A is given in (|l9|). It is easy to see that 

^;||C|| o<(-2|K||oo+7+llCl|oo)||CI|oo, 

so 1 1 C] |oo decreases when 

0< HClloo <2(l-M)- 7 , 
This implies that the basin of attraction of the sink contains this ball: 

D s {j) D {z : ||z — z*|| < 2(1 — M) — 7} 
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If the AI state s is to be in this basin, then we require that s £ D s (j), but since we know 
that ||z* — s||oo < Mqo, this is true when 

3M 00 ( 7 ) < 2 - 7 => 7 < 0.55566 . 

□ 



4 Conclusions 

A large class of dynamical systems can be formulated to have a nondegenerate AI limit. 
In this case a simple contraction mapping argument can be used to show that many of the 
symbolic orbits at the AI limit persist away from the limit. We have applied these ideas to 
the Henon map to show that all possible bounded orbits exist in the range (|7|). Translating 
this back to the original parameters k and b of the map (Q) gives 

k > 5 -±^(i + \b\r 



which is exactly the same bound as that found by Devaney and Nitecki [19| using a geomet- 
rical argument for the existence of a Smale horseshoe. The advantage of the AI argument 
is that it can be easily generalized to systems where the geometrical argument might be 
difficult, such as higher dimensional maps. Moreover, the AI theory can also be used to 



give bounds for the existence of orbits corresponding to various subshifts of finite type [23]. 

Our main interest in the AI limit is to use it as the basis for a numerical method to 
find orbits of various types. One such method we discussed here is the "pseudo-gradient" 
algorithm of Biham and Wenzel. We showed that this method, which was previously only 
justified heuristically, has a rigorous foundation close enough to an AI limit. Indeed, the 
signs in the matrix Diag(s) are determined by signature of the critical point at the AI limit. 
We also were able to show that the AI state can be used as an initial condition for the 
method. However, our theorem applies only for a limited range of parameter values, and 
one might reasonably apply the method over a wider range of parameters. 

Some care is advisable however, since it is known that the method of Biham and Wenzel 
can fail far from the AI limit We will show in a forthcoming paper [^] that in some 
applications, such as that in ||, the method does appear to find all orbits. 

Rather than use the Biham and Wenzel method we implement a simple "continuation" 



technique for our numerical studies [23|. This method is guaranteed to find all bounded 
orbits of the map that extend to the AI limit (there could be "bubbles" that do not extend, 
but as far as we know these have never been seen for the Henon case). For example, in the 
area preserving case (b = 1), our numerical results indicate that the horseshoe is destroyed 
at e = 0.41887923, when a pair of orbits homoclinic to the fixed point, s = {. . . , 1, 1, . . . } = 
{+°°}, collide in a saddle-node bifurcation. These orbits have the symbol sequences 

|_|_oo h°°} '■ > {+°° I h°°} 

This parameter value also corresponds to the accumulation point of an infinite number of 
periodic saddle-node bifurcations, and hence the first value at which the bounded orbits of 
the Henon map have nonzero measure. 
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